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ABSTRACT 



The astrometric data on the runaway star BD+43° 3654 are consistent with the origin of this 04If star in the center of the Cyg 0B2 
association, while BD+43° 3654 is younger than the association. To reconcile this discrepancy, we suggest that BD+43° 3654 is a 
blue straggler formed via a close encounter between two tight massive binaries in the core of Cyg 0B2. A possible implication of this 
suggestion is that the very massive (and therefore apparently very young) stars in Cyg OB2 could be blue stragglers as well. We also 
suggest that the binary-binary encounter producing BD-l-43° 3654 might be responsible for ejection of two high-velocity stars (the 
stripped helium cores of massive stars) - the progenitors of the pulsars B2020-I-28 and B2021+51. 
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1. Introduction 

BD-i-43° 3654 is a massive (04If) runaway star located ~ 2°7 
northeast of the Cyg OB2 association (Comeron & PasquaH 
12007 V The astrometric data on BD-i-43° 3654 suggest that this 
star was ejected from the center of CygOB2 ~ 1.8 Myr ago 
(Comeron & Pasquah |2007] |. BD-h43° 3654 is one of the three 
known very massive (^ 60 Mq) runaway stars (the two other 
stars are ^ Pup and A Cep) whose high (<; 40kms"') pecu- 
Uar velocities cannot be explained within the framework of the 
binary-supernova scenario (see Sect. 3). The more likely channel 
for production of these high-velocity massive objects is through 
dynamical processes in the dense cores of young massive star 
clusters (Poveda et al. 1967; Leonard & Duncan 1990), partic- 
ularly through close encounters between tight massive binary 
stars (Mikkola |1983l Leonard & Duncan 1990) . The most com- 
mon outcome of binary -binary encounters is the exchange of the 
more massive binary components into a new eccentric binary or 
a single merged star and ejection of the less massive ones with 
high velocities (e.g. Leonard [l995l l. The ejection velocity of the 
lightest components could be as high as the escape velocity from 
the surface of the most massive star in the binaries (i.e. ~ 1000 
kms"'; Leonard 1991). Correspondingly, the recoil velocity of 
the newly formed object (the massive binary or the merged star) 
could significantly exceed the escape velocity from the potential 
well of the parent cluster, so that this object becomes a runaway. 

In this Letter we suggest that BD-)-43° 3654 is a merged star 
formed via a close encounter between two tight massive bina- 
ries in the core of Cyg OB2 (cf. Leonard I1995I I. We suggest 
that BD-i-43° 3654 is a blue straggler, i.e. a rejuvenated star with 
the apparent age smaller than the age of the parent association 
(see Sect. 2 and Sect. 4). The high linear momentum imparted to 
BD-i-43° 3654 implies that the binary-binary encounter was very 
energetic and that the lower mass components of the binaries 
involved in the encounter were ejected with high velocities. 
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Interestingly, Cyg OB2 may be associated with the origin of 
two high-velocity pulsars, PSR B2020+28 and PSR B2021-H51. 
The high-precision proper motion and parallax measurements 
for these pulsars (presently separated by ~ 23°) lead to the sug- 
gestion that they originated ~ 1 .9 Myr ago within several par- 
sees of each other in the direction of Cyg OB2 (Vlemmings 
et al. 120041 ). Vlemmings et al. ( 120041 ) believe that the progen- 
itors of both pulsars were members of a common binary and 
that the pulsars were separated at the birth of the second one, 
following asymmetric supernova explosion. An alternative pos- 
sibility is that the pulsars were separated before their birth, 
and that they are the remnants of runaway stars ejected from 
CygOB2 due to the dynamical three- or four-body encounters 
(Gvaramadze 2007. Paper I). The relative position of the pulsars 
and BD-h43° 3654 on the sky and the similarity between their 
kinematic ages suggest that the three objects might be ejected 
from CygOB2 via the same dynamical event - the close en- 
counter between two massive binaries (Sect. 4; cf. Gies & Bolton 
1986i: Hoogerwerf et al. ,2001. Gualandris et al. |2004] l. 

2. BD+43° 3654 and Cyg 0B2 

The isochronal age of BD+43° 3654, derived by Comeron & 
Pasquali ^2007^ from the position of this 04If star in the 
Hertzsprung-Russell (HR) diagram and the isochrones from the 
evolutionary models for massive stars by Meynet et al. (1994), 
is similar to the kinematic age (i.e. the time since the ejection 
of the star from Cyg OB 2). Both ages should be compared with 
the age of the parent association of ~ 1-5 Myr (e.g. Bochkarev 
& Sitnik 1985; Torres-Dodgen et al. fTWT] Herrero et al. 19991 
Hanson 2003, Knodlseder |2003l l. The wide age spread could be 
considered as an indication that the star formation in Cyg OB2 
is non-coeval (Massey & Thompson 1 199 II Hanson l2003l l. On 
the other hand, the youngest ages of 5 3 Myr come from the 
presence in Cyg OB2 of a two dozen early type O stars, while 
the evolutionary status of the less massive stars in the associ- 
ation is consistent with an age of 4-5 Myr (the situation typi- 
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cal of young massive star clusters and associations; e.g. Massey 
12003 V One can hypothesize that the actual age of Cyg OB2 is ~ 
4-5 Myr and that the most massive (and therefore apparently the 
youngest) stars in the association are in fact the rejuvenated stars 
(blue stragglers) formed via merging of less massive stars in the 
course of close binary-binary encounters ~ 2 Myr ago, when 
CygOB2 was much more compact (see below). In our reason- 
ing we proceed from the results of binary-binary scattering ex- 
periments by Leonard ( |1995l l. which showed that about half of 
the merged stars formed via close encounters between binaries 
remain bound to the parent cluster (some of them form binaries 
with another blue straggler) and appear much younger than other 
members of the cluster (see also Portegies Zwart et al. |1999l ). 

An indirect support for our hypothesis comes from the re- 
cent study of early-A stars in the direction of Cyg OB2 by Drew 
et al. ( I2OO8I 1. This study revealed several hundred A stars within 
the boundaries of the association and suggested that the age of 
these stars is > 5 Myr, provided that they are located at a dis- 
tance of < 1.7 kpc (i.e. at the distance of CygOB2; Kiminki et 
al. 120071 cf. Torres-Dodgen et al. 1991; Massey & Thompson 
II99II 1. The age of ~ 3-5 Myr would also be required if: (i) the 
pulsars B2020-h28 and B2021+51 indeed originate in CygOB2 
(Paper I; see also Sect. 4), (ii) the HEGRA TeV source observed 
in the direction of Cyg OB2 (Aharonian et al. l2002l l is related to 
a young pulsar in the association (Bednarek |2003l l and (iii) a par- 
tially non-thermal shell-like object coincident with the HEGRA 
source is a supernova remnant (Butt et al. l2008l l. 

Moreover, the stellar content of Cyg OB2 could be contami- 
nated by young massive stars injected into the association from 
nearby associations (Uyaniker et al. 200l]l and the numerous 
young clusters around Cyg OB^] (Dutra & Bica l200lt Comeron 
& Torra |200Tt Le Duigou & Knodlseder |2002l ). In the latter case, 
some of the 'alien' O stars can be produced through merging 
of two or three stars in the course of binary-binary encounters 
in clusters of B-type stars (see Leonard 1995). For example, 
in our search for bow shocks around Cyg 0B2 (Gvaramadze & 
Bomans, in preparation), we discovered a bow shock produced 
by one of the early-type stars from the list of new members of the 
association by Comeron et al. (2002). This star [designated by 
Comeron et al. ( 120021 1 as A37] was classified by Hanson (I2003I I 
as a 05V star, and therefore it should be a young (^ 1 Myr) 
object. The photometric distance to this star of 1.7-1.8 kpc [de- 
rived with use of the UBVJHK synthetic photometry of Galactic 
stars by Martins & Plez (20()6.)] is consistent with the distance 
to CygOB2. The astrometric data on the star and the geome- 
try of the bow shock, however, suggest that this high-velocity 
(~ 120 km s"') runaway was ejected either from the young clus- 
ter embedded in the H II region DR 15 (located on the border of 
CygOB2 ~ 1° west of the current position of the star) or from 
the open cluster NGC 6913 (centered ~ 3°A west of the star). 

At the distance to Cyg OB2 of ~ 1 .7 kpc, the peculiar (trans- 
verse) velocity of BD-i-43° 3654 is ^ 40 + lOkms"' [we used 
here the Galactic constants Rq - kpc and ©o = 200 km s ' 
(e.g. Reid 1993; Kalirai et al. |2004l Avedisova |2005] l and the 
solar peculiar moti on (U q, Vq, Wq) = (10.00,5.25,7. 17) kms"' 
(Dehnen & BinneyTg^F); cf. Comeron & Pasquali ("2007)]. The 
position of BDh-43° 3654 in the HR diagram and the evolution- 
ary tracks by Meynet et al. ( |1994| l imply an initial mass of the 



' A study of the Cygnus X region by Schneider et al. ( 12006) suggests 
that these clusters and the Cyg 0B2 and 0B9 associations form a co- 
herent complex and that formation of the clusters was triggered by the 
effect of massive stars in the associations (cf. Le Duigou & Knodlseder 
|2002t . 



star of ^ 70 + 15 Mq (Comeron & PasquaU |2007] l. The high Un- 
ear momentum attained by this runaway star could be used to 
constrain the possible mechanisms of its origin. In Sect. 3 we 
show that runaways of this mass and velocity are unlikely to be 
produced via the disruption of a binary due to the (asymmetric) 
supernova explosion (the binary-supernova scenario; Blaauw 
[1961; Stone 1991). Another possibility is that this star attained 
a high peculiar velocity via the strong dynamical three- or four- 
body encounter (the dynamical ejection scenario; Poveda et al. 
1967, Gies & Bolton 1986). In Sect. 4, we suggest that the most 
likely path for the origin of BD-i-43° 3654 is through the close 
encounter between two tight massive binaries. 

Cyg OB2 is one of the most compact and massive associa- 
tions in the Milky Way. It contains ~ 100 O stars or stars with 
O-type progenitors (Knodlseder 2000, Comeron et al. I2002I I. 
The half Ught radius of CygOB2 is ~ 6 pc (Knodlseder |2000] l. 
Assuming that the association expands with a velocity equal to 
its velocity dispersion (~ 2.4 km s"'; Kiminki et al. 2007), one 
finds that the majority of massive stars in Cyg OB2 were origi- 
nally concentrated in a region of radius of < 1 pc (that is consis- 
tent with the observation that the initial radii of young clusters 
are ^ 1 pc; Kroupa & Boily 2002 ). It is therefore plausible that at 
the moment of ejection of BD-i-43° 3654 the stellar number den- 
sity in the core of Cyg 0B2 was high enough to ensure that close 
encounters between its constituents were frequent. The neces- 
sary condition for effective production of runaways is a high 
binary fraction among massive stars. The recent radial veloc- 
ity survey of CygOB2 by Kiminki et al. ( 120071 ) gives a lower 
limit on the massive binary fraction of 30 - 40%, while the com- 
parison of the data from the survey with the expectations of the 
Monte Carlo models suggest that this fraction could be > 80% 
(Kobulnicky & Frver l20(T7l l. Thus we believe that ~ 2 Myr ago 
the conditions in the core of Cyg OB2 were favourable for the 
dynamical processes discussed in Sect. 4. 



3. BD+43° 3654: binary-supernova scenario 

According to the binary-supernova scenario, a massive star re- 
siding in a binary system could attain a high peculiar velocity 
due to the disruption of the binary after the companion star ex- 
ploded as a supernova (Blaauw 119611 Stone |1991| l. We show 
that the peculiar velocity of BD+43° 3654 cannot be accounted 
for within the framework of this scenario (cf. Vanbeveren et al. 
I2007T I. 

The kinematic age of BD-h43° 3654 of ~ 1.8 Myr and the 
minimum possible lifetime of the supernova progenitor star of ~ 
2.5-3 Myr imply that the actual age of BD-i-43° 3654 should be ^ 
4.5-5 Myr. The discrepancy between the 'observed' and the in- 
ferred ages could be reconciled if BD+43° 3654 is a rejuvenated 
star. For example, it could be rejuvenated through mass transfer 
from the primary star during the Roche lobe overflow stage (e.g. 
Dray & Tout 2007 ). In this case, the effect of rejuvenation would 
be significant only if the mass gained by the rejuvenated star 
was larger than its initial mass. Although one cannot exclude 
this possibility, we note that stars more massive than ~ 40 M© 
could lose a significant fraction of their mass via the heavy stellar 
wind and the Roche lobe overflow will not occur (Vanbeveren et 
al. ,1998i ). Another possibiHty is that BD+43° 3654 was formed 
in the course of a close encounter between two massive bina- 
ries, during which two stars of the binaries merged into a single 
rejuvenated star (now BD+43° 3654) and caught a third (more 
massive or more evolved) star to form a new binary, while the 
fourth star was ejected as single (cf. Leonard|T995]l. 
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Let us assume that at the moment of supernova explosion in a 
binary system the mass of the second star was 70 M©. It is obvi- 
ous that to disrupt such a massive system, the explosion should 
be asymmetric so that the supernova stellar remnant attained a 
kick. In this case, the stellar remnant can impart some momen- 
tum to the companion star in the course of disintegration of the 
binary (Tauris & Takens 1998'). The magnitude of the momen- 
tum depends on the angle between the kick vector and the direc- 
tion of motion of the exploding star and reaches its maximum 
for a certain value of the angle [given by Eq. (4) in Gvaramadze 
120061 . One can show that to accelerate a star as massive as 
BD-i-43°3654 to the velocity of 40kms"', the kick direction 
should be very carefully tuned (i.e. should be within several de- 
grees of the direction towards the companion star), while the 
magnitude of the kick should be very large (> 700kms ', if the 
stellar remnant is a neutron star, or > 200 km s"', if the rem- 
nant is a black hole of mass ~ 5 Mq). Although one cannot 
exclude this possibility, we consider it as highly unlikely (cf. 
Gvaramadze 2006; Paper I). 
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Fig. 1. Relative position of BD+43°3654 and the pulsars 
PSR B2020-H28 and PSR B2021-H51 on the sky. The circle of 
angular diameter of 2° indicates the boundaries of the Cyg OB2 
association. 



4. BD+43° 3654: dynamical ejection scenario 

The numerical experiments by Leonard ( 1 19951 1 showed that 
a significant fraction of unbound blue stragglers (formed via 
binary-binary encounters in dense clusters) attain peculiar veloc- 
ities large enough (> 30kms"') to be classified as runaways. It is 
therefore tempting to consider the possibility that BD-i-43° 3654 
is a blue straggler formed through an encounter between two 
tight massive binaries. The most common outcome of encoun- 
ters between such binaries is the exchange of the more massive 
components into a new eccentric binary or a single merged star 
and ejection of the less massive ones with high velocities. The 
binary would ultimately coalesce into a single star if its orbit is 
sufficiently compact. 

Let us assume that BDh-43° 3654 is the result of a merging of 
two main-sequence (MS) stars of mass M\ and M2 {Mi > M2) 
and of the same age f^erg - 3 Myr. Under these assumptions, 
the merger product is also an MS star (of mass M» ^ Mi -1- M2), 
and its new age is given by (see Meurs & van den Heuvel [T989l 
Portegies Zwart et al. |1999l ): 

Ml tusiM.) ^ 
'*~M. fMs(MO 

where tus (M.) and tusiMi) are the MS lifetimes of stars of 
mass M» and Mi . One can consider two cases: (i) BD-h43° 3654 
is the product of a physical collision of two stars during the close 
binary-binary encounter and (ii) BD-i-43° 3654 is the result of co- 
alescence of two stars in a close binary system. In the first case, 
one should require that f, ~ 0, since the current (apparent) age 
of BD-i-43° 3654, t - + fidn, is comparable to its kinematic 
age fkin - 1.8 Myr This requirement could be fulfilled only 
if the encounter between binaries occurs very soon after their 
birth in the association (i.e. f,nerg ~ Myr), which contradicts 
our assumption that f^erg ~ 3 Myr. Note, however, that here we 
neglected the possibility that the colliding stars were already re- 
juvenated through mass transfer in the original tight binaries so 
that their apparent age was much less than the actual one. In the 
second case, one can assume that the binary components coa- 
lesced into a single star only recently, i.e. the time elapsed since 
the formation of BD+43° 3654 is less than the kinematic age 
of the binary. As an example, suppose that the runaway binary 
consists of two stars of mass ~ 35 Mq. For Jms (35 M©) - 4.5 
Myr and fMS (70 Mq) ^ 2.5 Myr (Meynet & Maeder |20U3l l, and 



assuming that the binary merged ~ 1 Myr after ejection from 
Cyg OB2 (i.e. fmeig ~ 4 Myr), one has from Eq. (1) that ~ 1.1 
Myr During the next ~ 0.7 Myr (= fyn - ?*), the merged star 
evolves into a blue supergiant with parameters similar to those 
of BD-i-43° 3654 [see the evolutionary tracks by Meynet et al. 
( 119941 ) and the calibration of parameters of Galactic O stars by 
Martins et al. (|2005l l1. 

In Sect. 1, we mentioned that the origin of two pulsars, 
B2020+28 and B2021-H51, could be associated with CygOB2. 
In Paper I, we suggested that these pulsars are the remnants of 
runaway stars ejected (with velocities similar to those of the pul- 
sars) from Cyg OB2 due to the dynamical three- or four-body en- 
counters. Our suggestion was based on the recently recognized 
fact that the high-velocity pulsars could be the descendants of 
high-velocity runaway stars (i.e. the peculiar velocities of pulsars 
do not necessarily originate from asymmetric supernova explo- 
sions; Paper I; Gvaramadze et al. 120081 1. Strong support for this 
possibility comes from the discovery of early B-type stars mov- 
ing with velocities of ~ 500 -700 km s"' (Edelmann et al. 120051 
Przybilla et al. |2008l Heber et al. '2008^. In Paper I, we used 
the similarity between the spin-down and the kinematic ages of 
B2020-I-28 and B2021+51 to suggest that their progenitors were 
the short-lived (< 1 Myr) helium cores of massive stars, while 
from the age of Cyg OB2 at the moment of ejection of the he- 
lium cores (~ 3 Myr) we inferred that the zero-age MS masses 
of the ejected stars were > 50-60 M©; it is believed that stars of 
this initial mass could leave behind a neutron star (see Woosley 
et al. 1995; Muno et al.'2006', Bibby et al. |2008] l. The relative 
position of the pulsars and BD+43° 3654 on the sky (see Fig. 1) 
and the similarity between their kinematic ages suggest that the 
three objects might have had a common origin in the close en- 
counter between two massive binaries. Below, we discuss this 
possibility in detail. 

The (transverse) recoil velocity attained by the merged star 
is given by: 

V.,tr = M;' [(mi Vi^af + (OT2V2,tr)^ + 2mim2Vl,trV2,tr COS Or] ' , (2) 

where mi, mi and vi ^r, V2,ti are, respectively, the masses and the 
transverse velocities of the helium cores (the progenitors of the 
pulsars) and a is the angle between the ejection velocity vectors 
of the helium cores. We assume that the pulsars did not received 
(significant) kicks at birth and therefore move (almost) in the 
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same direction as their progenitors. Adopting that the binary- 
binary encounter occured 1.8 Myr ago and neglecting the effect 
of the Galactic gravitational potential, one obtains the transverse 
velocities of B2020+28 and B2021+51 as vi,tr - 200kms-' 
and V2,tr - 170 km s"', and a - 160°. Using Eq. (2), one 
can show that to explain the 'observed' transverse velocity of 
BD+43°3654 of 40kms-', the progenitor of B202 1+51 should 
already lose most of its mass via the stellar wind, i.e. by the 
moment of binary-binary collision 1112 should be several times 
smaller than mi. For M, = 70 + 15Mq and assuming that mi - 
(3 - 4)ot2, one has from Eq. (2) that mi ^ (17 - 19) + 4 Mq (i.e. a 
quite reasonable figure). For the above parameters, one finds that 
the velocity vector of the recoiled merged star is somewhat mis- 
ahgned (~ 10°) to the residual velocity vector of BD+43° 3654. 
To explain this misalignment, one can assume that Cyg OB2 has 
a peculiar (transverse) velocity of ~ 7kms"' in the northwest 
direction (cf. Hoogerwerf et al. 120011) . The origin of peculiar ve- 
locity of this magnitude (typical of the OB associations near the 
Sun; de Zeeuw et al. 1 19991 1 could be understood if formation of 
Cyg OB2 was triggered by the collision between two molecular 
clouds (cf. Schneider et al. 2006). 

One can also constrain the radial velocity of BD+43° 3654 
using the parallactic distances to B2020+28 and B2021-H51, re- 
spectively, of 2.7+1;] ^-^-02 (Vlemmings et al. |20D4l ). 
Taken at face value, these distances imply the pulsar radial ve- 
locities vi,r ^ 530!®°kms-i and V2,,- ^ 160![^°kms-', while 
from the conservation of the linear momentum one has V* ^ ^ 
-(150^jjg)kms"'. The peculiar radial velocity of CygOB2 
of ^ -6kms"' (derived from the mean systemic velocity of 
the association of ^ -lOkms"'; Kiminki et al. 2007) only 
slightly changes the above figures. Our scenario for the origin 
of BD+43° 3654 therefore suggests that the radial component of 
the peculiar velocity of this star should be negative and at least 
as large as the transverse one. It could be that the radial velocity 
measurements for BD+43° 3654 will invalidate our scenario, so 
that the origin of the runaway star and the pulsars would not be 
related to each other Even in this case, we believe that the origin 
of BD-h43° 3654 should be accompanied by the ejection of two 
high-velocity stars, either early type B-stars or stripped helium 
cores of the more massive stars. 

Thus, we suggest that the runaway massive star 
BD-i-43° 3654 originate from a close encounter between 
two binaries, originally consisting of a ~ 35 and a ~ 50 - 60 Mq 
star. The more massive stars in each system evolved for ~ 3 
Myr, losing most of their mass, becoming helium stars. At this 
point, the binaries interacted, ejecting the helium stars at high 
velocity and resulting in the merger or near-merger of the 35 M© 
stars into a 70 M© star, which recoiled at a proportionally lower 
velocity. The more massive object is now seen as a blue straggler 
04If star, while the helium stars exploded as supernovae soon 
after the ejection and produced the pulsars with a small or no 
kick at birth. 
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